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End  plates  of  the  proton  exchange  membrane  fuel  cell  (PEMFC)  need  to  be  well  designed  because  their 
strength  and  rigidity  directly  affect  the  clamping  pressure  distribution  and  thus  affect  the  performance 
and  lifetime  of  fuel  cell  stacks.  In  this  paper,  a  multi-objective  topology  optimization  model  of  the  end 
plates  in  a  PEMFC  stack  with  nonlinear  contact  boundary  conditions  was  established  to  obtain  an  opti¬ 
mized  structural  design.  It  was  found  that  the  design  improved  with  topology  optimization  is  not  only 
light  but  also  meets  manufacturability  requirements.  This  provides  good  guidance  for  the  design  of  a 
high-performance  end  plate. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  has  attracted 
much  attention  as  a  promising  power  source  with  wide  applica¬ 
tion  owing  to  its  high  power  efficiency,  high  energy  density,  quick 
cold-start  capability  and  low  level  of  pollution  [1-4].  It  is  a  possible 
substitute  for  combustion  engines  in  automobiles  and  heat  gener¬ 
ators  in  houses.  For  PEMFCs  to  meet  the  requirements  of  practical 
applications,  their  performance  and  lifetime  need  to  be  optimized. 

As  the  voltage  and  total  power  generated  by  a  single  PEMFC 
is  rather  limited,  PEMFC  products  for  practical  applications  com¬ 
prise  numerous  planar  single  cells  in  series  that  form  large  stack 
structures  with  a  cyclic  cell  characteristic.  Typically,  two  end  plates, 
between  which  the  unit  cells  are  inserted,  are  located  on  the  out¬ 
side  of  the  stack  to  provide  proper  internal  compression,  and  the 
entire  fuel  cell  structure  is  fastened  with  a  set  of  bolts  or  by  another 
fastening  mechanism. 

End  plates  need  to  be  well  designed  to  have  sufficient  mechani¬ 
cal  strength  and  stiffness  while  consisting  of  a  minimum  amount  of 
materials.  A  good  design  of  the  PEMFC  stack  requires  that  the  con¬ 
tact  pressure  in  each  cell  is  not  only  approximately  equal  but  also 
has  as  uniform  distribution  as  possible  [5].  In  general,  there  may 
be  high  contact  electrical  resistance  at  the  interface  between  the 
gas  diffusion  layer  (GDL)  and  the  bipolar  plate  (BPP)  if  the  contact 
pressure  is  unreasonably  small  [6-8].  However,  the  permeability 
of  the  GDL  will  be  too  low  to  provide  high  reaction  efficiency  if  the 
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contact  pressure  is  excessively  large  [9].  In  the  practical  design  of 
a  fuel  cell  stack  system,  the  end  plates  are  designed  to  be  as  rigid 
as  possible  so  that  the  clamping  load  applied  by  the  bolt  assembly 
can  be  uniformly  distributed  within  the  whole  fuel  cell  stack. 

From  a  design  point  of  view,  high  rigidity  of  the  end  plate 
can  be  simply  achieved  by  increasing  its  thickness.  However,  an 
end  plate  that  is  too  thick  will  provide  a  too  large  mass  or  vol¬ 
ume  [5].  To  achieve  low  mass  and  high  rigidity  simultaneously, 
end  plates  with  complex  configurations,  such  as  ribbed-,  bomb- 
,  or  bow-shaped  end  plates  [10],  have  been  developed.  Although 
most  of  these  designs  seem  to  be  fairly  effective,  none  are  based  on 
topology  optimization,  a  powerful  theoretical  tool  for  designers  to 
find  the  best  topology  or  layout  for  given  objectives  under  specific 
constraints  [11]. 

In  this  paper,  a  multi-objective  topology  optimization  model 
based  on  finite  element  analysis  (FEA)  is  proposed  to  obtain  the 
optimized  material  distribution  (topology)  of  the  end  plate  in  a  spe¬ 
cific  design  space.  The  optimization  objectives  are  to  maximize  the 
stiffness  of  the  end  plate  and  unify  the  pressure  distribution  at  the 
contact  interfaces  of  the  fuel  cell  stack  at  the  same  time.  Combined 
with  the  equivalent  stiffness  model  proposed  in  our  previous  works 
[  1 2,1 3  ],  the  optimization  design  is  realized  and  verified  respectively 
using  the  commercial  software  Optistruct®  and  Abaqus®. 

2.  Basic  theory 

2.2.  Topology  optimization 

According  to  the  types  of  design  variables,  the  optimization 
design  of  structures  (only  continuum  structures  are  considered  in 
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Fig.  1.  Schematic  diagram  of  the  relationship  between  the  optimization  methods 
and  design  stages. 


this  paper)  can  be  divided  into  three  levels  -  size,  shape  and  topol¬ 
ogy  optimizations  -  corresponding  to  the  three  stages  of  the  design 
process  [11],  as  shown  in  Fig.  1.  Topology  optimization  is  carried 
out  at  the  top  level  and  aims  to  find  the  best  distribution  of  mate¬ 
rial  in  the  specified  design  space,  thus  setting  the  optimized  route 
for  load  transmittance.  The  importance  of  this  optimization  lies  in 
the  fact  that  the  choice  of  the  appropriate  topology  of  a  structure 
in  the  conceptual  design  stage  is  generally  the  most  decisive  factor 
of  the  product’s  efficiency.  Moreover,  size  and  shape  optimizations 
do  not  change  the  structural  topology  in  the  process  of  finding  the 
optimal  solution.  Therefore,  the  topology  optimization  is  a  valuable 
tool  in  preprocessing  of  the  size  and  shape  optimizations  [14,15]. 

Since  the  homogenous  method  was  first  introduced  by  Bendsoe 
and  Kikuchi  [16]  in  1988  to  solve  topology  optimization  problems 
of  continuum  structures,  many  new  methods  have  been  developed, 
such  as  the  variable  density  method  (VDM),  variable  thickness 
method,  and  evolutionary  structural  optimization.  Among  the 
methods,  the  VDM  is  one  of  the  most  widely  used  and  is  often  inte¬ 
grated  with  FEA.  In  the  VDM,  the  material  density  of  each  element  is 
considered  as  the  design  variable  and  varies  continuously  between 
0  and  1 ,  where  “0”  represents  a  “void”  and  “1”  represents  the  “solid”. 
Intermediate  values  of  the  density  between  0  and  1  represent  fic¬ 
titious  material.  Originally,  the  stiffness  of  material  is  assumed  to 
depend  linearly  on  the  density.  However,  this  assumption  usually 
causes  numerous  gray  areas  in  the  intermediate  density  region, 
which  are  impracticable  when  determining  the  topology  of  a  given 
material  or  using  different  materials  within  the  design  space. 
Hence,  a  technique  is  needed  to  penalize  the  intermediate  den¬ 
sities  and  force  the  final  design  density  to  be  approximately  either 
0  or  1.  In  practice,  the  solid  isotropic  microstructure  with  penalty 
(SIMP)  method  [17-19]  is  widely  used.  It  can  be  expressed  for  any 
two-dimensional  or  three-dimensional  solid  element  as 


K'{p)  =  ppxK,  (1) 

where  I<'  and  I<  are  respectively  the  penalized  and  real  stiffness 
matrixes  of  an  element,  p  is  the  relative  density  (a  design  vari¬ 
able),  and  P  is  the  penalization  factor.  The  relationship  among  these 
parameters  is  shown  in  Fig.  2.  It  is  found  that  the  larger  the  penal¬ 
ization  factor,  the  more  discrete  the  obtained  results.  Bendsoe  and 
Sigmund  [20]  showed  that  good  results  can  be  obtained  using  the 
SIMP  method  when  P  >  3.  In  this  paper,  the  SIMP  method  with  P= 3 
is  employed  for  local  topology  optimization  problems. 


2.2.  Multi-objective  optimization 

In  a  practical  optimization  problem,  there  are  usually  several 
conflicting  design  criteria.  For  example,  objectives  to  reduce  the 
weight  and  increase  the  stiffness  of  a  structure  cannot  simply  coex¬ 
ist  in  one  optimization  problem  since  only  a  single  objective  can 
be  considered  in  each  optimization  process.  Consequently,  a  so- 
called  multi-objective  (multi-criteria,  vector,  Pareto)  optimization 
method  has  been  developed  as  a  methodology  for  solving  opti¬ 
mization  problems  with  several  simultaneous  objective  functions 
[21]. 

In  general,  multiple  objectives  need  to  be  transformed  into  a 
single  objective  by  retaining  one  selected  criterion  to  be  optimized 
while  treating  the  remaining  objectives  as  either  constraints  or  by 
assigning  suitable  weights.  However,  the  choice  of  constraint  limits 
may  be  a  difficult  task  because  most  design  quantities  do  not  have  a 
constraint  nature.  Thus,  a  norm  method  with  weighting  coefficients 
[21-23]  is  proposed  in  this  paper. 

2.3.  Equivalent  stiffness  model  of  the  PEMFC  stack 

An  equivalent  stiffness  model  of  the  PEMFC  stack  has  been 
established  by  the  authors  [12,13].  Using  this  model,  the  whole  fuel 
cell  stack  is  simplified  as  a  mechanical  model  consisting  of  a  large 
number  of  elastic  elements  (springs)  in  either  series  or  parallel. 
Thus,  the  stack  structures  in  this  paper,  except  for  the  end  plates, 
insulator  plates,  and  current  collectors,  can  be  reduced  to  cuboids 
with  specific  stiffness  to  reduce  the  task  of  topology  optimization. 
Young’s  modulus  of  the  cuboids  can  be  calculated  using 


where  k  is  the  equivalent  stiffness  of  the  corresponding  component 
[12],  E  is  Young’s  modulus,  A  is  the  cross-section  area,  and  L  is  the 
length  along  the  axial  direction. 

3.  Model  demonstrations 

3.1.  Geometry 

A  three-dimensional  FEA  model  including  the  insulator  plates 
and  current  collectors  is  established  on  the  basis  of  a  previ¬ 
ously  designed  PEMFC  stack  structure  [24]  but  with  the  following 
changes  and  simplifications: 

(a)  Parallel  flow  channel  structure  is  one  of  the  most  com¬ 
mon  structures  used  in  PEMFC  stacks  owing  to  its  simplicity 
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Fig.  3.  Dimensions  of  the  single-cell  structure. 


in  design  and  processing.  Therefore,  it  is  adopted  in  this 
paper. 

(b)  Because  thermal  effect  is  not  taken  into  account  in  this  paper, 
existence  of  the  cooling  channels  might  only  affect  the  struc¬ 
tural  stiffness.  However,  according  to  the  equivalent  stiffness 
model  [12,13],  the  cooling  channel  structure  with  a  small 
cross-section  has  a  minor  influence  on  the  structural  stiff¬ 
ness.  Although  any  types  of  the  cooling  channel  structure  can 
be  included  in  our  optimization  model,  the  original  cooling 
channels  in  BPPs  will  not  considered  in  the  present  paper  for 
simplicity. 

(c)  Gasket  sealants  like  ones  given  in  Ref.  [25]  are  adopted  because 
they  show  roughly  a  linear  mechanical  response  under  clamp¬ 
ing  load  and  have  been  widely  used  in  practical  applications. 

In  the  following,  the  single  PEMFC  structure  with  six-bolt  assem¬ 
bly  similar  to  that  in  the  previous  design  [24],  as  shown  in  Fig.  3,  is 
studied.  Only  1  /8  of  the  single  cell  is  considered  owing  to  structural 
symmetry.  The  thicknesses  of  the  PEM  and  GDL  are  respectively 
0.05  and  0.275  mm  [25-27]. 

3.2.  Finite  element  mesh  and  material  properties 

A  mapped  meshing  method  is  adopted  to  ensure  proper  element 
connectivity  and  a  reasonable  aspect  ratio.  Moreover,  there  are 


restrictions  to  the  element  size  in  the  design  area.  Since  extremely 
small  features  produced  in  the  topology  optimization  are  usually 
unacceptable  in  the  manufacturing  process,  a  minimum  dimen¬ 
sion  constraint  should  be  introduced  in  the  optimization  model. 
Usually,  the  mean  element  size  should  be  less  than  one-third  of 
the  minimum  dimension  constraint.  As  the  end  plates  studied  here 
are  made  of  aluminum  alloy,  the  minimum  dimension  constraint 
is  assumed  to  be  6  mm  according  to  the  corresponding  machin¬ 
ing  process.  Therefore,  the  mean  element  size  for  the  end  plates  is 
chosen  to  be  2  mm. 

Material  properties  of  each  component  in  the  stack  are  listed  in 
Table  1 .  All  components  are  assumed  to  have  linear  elastic  behavior. 

3.3.  Boundary  conditions  and  optimization  parameter  settings 

Symmetry  constraints  are  applied  to  the  symmetrical  planes  to 
prevent  rigid-body  displacements.  A  contact  algorithm  is  used  to 
simulate  the  interactions  between  the  two  components  as  shown  in 
Fig.  4.  Therefore,  highly  nonlinear  factors  are  involved  in  this  topol¬ 
ogy  optimization  model.  The  clamping  load  is  calculated  from  the 
tightening  torque  (16  Nm)  [24]  and  applied  at  the  corresponding 
positions  of  the  bolts. 

The  design  area  should  be  specified  before  setting  the  optimiza¬ 
tion  model.  The  material  of  the  end  plate  far  from  the  clamping 
bolts  is  removed  as  it  contributes  little  to  the  assembly  load  bear- 


Table  1 

Material  properties  of  the  stack  hardware. 
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Component 

Material 

Young’s  modulus  (GPa) 

Poisson’s  ratio 

Source 

PEM 

Nafion® 

0.19 

0.25 

[26] 

GDL 

Carbon  paper 

10 

0.25 

[26] 

Bipolar  plate 

Graphite 

10 

0.25 

[26] 

Sealant 

VMQ. 

5.5 

0.3 

[28] 

Current  collector 

Copper 

100 

0.33 

[29] 

Insulation  plate 

POM 

2.6 

0.386 

[29] 

End  plate 

Aluminum  alloy 

69 

0.33 

[29] 

Clamping  bolt  and  nut 

SS  304 

209 

0.28 

[29] 

Washer 

Aluminum  bronze 

110 

0.3 

[29] 

contact 


equivalent  stiffness  block 


end  plate 
(design  area) 


end  plate 
(non-design  area) 

/ 


contact 


contact 


insulator  plate 


/ 

cum 

/ 


current  collector 


Fig.  4.  Schematic  diagram  of  the  optimization  model.  (Only  the  mesh  of  the  design 
area  is  displayed.) 


ing.  Since  the  topology  optimization  is  sensitive  to  the  initial  design 
area,  which  should  be  as  large  as  possible  [11],  only  regions  around 
the  bolt  holes,  the  supported  edge  and  the  bottom  side  are  excluded 
during  optimizing  (see  Fig.  4). 

There  are  two  objectives  that  need  to  be  considered  simulta¬ 
neously:  maximizing  the  structural  rigidity  of  the  end  plate  and 
uniformizing  the  pressure  distribution  in  the  whole  stack.  The 
former  can  be  carried  out  through  minimizing  the  structural  com¬ 
pliance  [23]  and  the  latter  can  be  achieved  by  minimizing  the  root 
mean  square  (RMS)  of  the  displacements  at  all  nodes  in  the  con¬ 
tact  area  between  the  current  collector  and  the  equivalent  stiffness 
block.  Consequently,  the  multi-objective  of  the  model  is  formulated 
as  [21-23] 


minimize  : 


co. 


C(p)  -  C, 


min 


c, 


L  '-max 


c, 


mm 


+  00- 


RMS{p)  -  RMS 


min 


RMS 


max 


RMS 


min 


i/q 


,(3) 


where  oo\  and  co2  are  weighting  coefficients,  C(p)  and  RMS{p) 
respectively  are  functions  of  the  compliance  of  the  end  plate  and 
the  root  mean  square  mentioned  above,  Cmax  and  RMSmax  are  the 
corresponding  responses  determined  for  the  original  structure,  and 
Cmin  and  RMSm[n  are  calculated  by  solving  the  single-objective  opti¬ 
mization  problems  expressed  as 


minimize  :  C(p),  (4) 

and 


minimize  :  RMS(p).  (5) 

Here,  oo\  =  oo2  =  0.5  and  q  =  0.5  are  chosen  for  this  problem.  The  con¬ 
straint  of  this  optimization  model  can  be  written  as 

Volume  fraction  =  — — ^non~design  <  0.3,  (6) 

''design 

where  V  is  the  total  volume  of  the  whole  model  in  the  current  iter¬ 
ation,  Vnon_design  is  the  initial  non-design  volume,  and  Vdesign  is  the 


initial  design  volume.  The  objective  is  compiled  in  FORTRAN  lan¬ 
guage  as  user-defined  external  functions.  The  model  built  here  is 
referred  to  as  the  reference  optimization  model  in  the  following. 

4.  Optimization  results 

Results  of  the  topology  optimization  are  usually  expressed 
by  the  relative  density  of  each  element.  Therefore,  the  elements 
with  low  densities  are  hidden  to  provide  clear  guidance  for  the 
establishment  of  the  optimized  topology  of  the  structure.  In  the 
following  figures,  the  color  difference  expresses  the  morphol¬ 
ogy  of  end  plate.  If  the  single-objective  optimization  method  is 
used,  to  increase  the  rigidity  of  the  end  plate,  material  should  be 
concentrated  in  the  area  linking  the  two  adjacent  bolt  holes,  as 
shown  in  Fig.  5(a).  However,  this  will  definitely  increase  elastic 
deformation  in  the  center  of  the  end  plate,  and  thus,  the  con¬ 
tact  pressure  in  the  GDLs  will  be  terribly  non-uniform.  On  the 
other  hand,  to  obtain  a  pressure  distribution  as  uniform  as  pos¬ 
sible  at  the  contact  interfaces,  it  is  effective  to  stiffen  the  structure 
in  the  contact  area  as  depicted  in  Fig.  5(b).  Yet  in  this  case,  there 
may  be  large  deformation  and  unreasonable  stress  concentration 


Fig.  5.  Schematic  diagram  of  the  single-objective  optimization  designs:  (a)  minimiz¬ 
ing  the  structural  compliance  and  (b)  minimizing  the  RMS  of  the  displacements  at  all 
nodes  in  the  corresponding  contact  area.  Elements  displayed  here  (and  hereinafter) 
have  relative  densities  greater  than  0.3. 
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Fig.  6.  Schematic  diagram  of  the  multi-objective  optimization  design:  (a)  original  end  plate;  (b)  optimization  result;  (c)  manufacturing  design  according  to  the  optimization 
result;  (d)  contact  pressure  distribution  on  the  GDL  for  the  original  structure;  (e)  contact  pressure  distribution  on  the  GDL  for  the  optimized  structure. 


in  the  end  plate,  which  are  obviously  against  the  design  pur¬ 
pose. 

Fig.  6  shows  the  detailed  multi-objective  optimization  process 
for  the  end  plate.  Fig.  6(a)  is  the  original  end  plate.  Fig.  6(b)  shows 
the  optimized  topology  considering  the  multi-objective,  which  can 
be  approximately  regarded  as  the  combination  of  results  depicted 
in  Fig.  5(a)  and  (b).  A  new  geometry  of  the  end  plate  is  obtained  (see 
Fig.  6(c))  according  to  the  material  distribution  given  in  Fig.  6(b). 
The  optimized  structure  is  validated  using  nonlinear  contact  FEA 
models  established  in  Abaqus®.  It  is  found  that  when  the  fuel  cell 


Fig.  7.  Schematic  diagram  of  the  multi-objective  optimization  designs  for  different 
numbers  of  clamping  bolts:  (a)  4  clamping  bolts  and  (b)  10  clamping  bolts. 


stack  is  equipped  with  the  original  end  plate,  the  contact  pressure 
at  the  center  of  the  GDL  is  much  less  than  that  at  the  edge,  as  shown 
in  Fig.  6(d).  Hence,  the  pressure  distribution  on  the  contact  surface 
between  the  GDL  and  BPP  is  obviously  non-uniform,  which  may 
induce  an  increase  in  the  interfacial  ohmic  resistance  in  the  cen¬ 
ter  region.  On  the  other  hand,  if  the  optimized  structure  is  used 
instead,  the  contact  pressure  distribution  over  the  whole  cell  area 
is  very  much  improved  (Pig.  6(e)).  Hence,  the  contact  pressure  with 
optimization  is  much  better  than  that  without.  Moreover,  the  vol¬ 
ume  of  the  optimized  end  plate  as  shown  in  Pig.  6(c)  is  reduced 


Fig.  8.  Schematic  diagram  of  the  multi-objective  optimization  designs  for  different 
numbers  of  cells  (six-bolt  assembly):  (a)  5-cell  stack  and  (b)  10-cell  stack. 
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Fig.  9.  Schematic  diagram  of  the  multi-objective  optimization  designs  for  several  combinations  of  weighting  coefficients:  (a)  co\  =0.7,  co2  =0.3;  (b)  co\  =0.6,  co2  =  0.4;  (c) 
( 0\  =  0.4,  o.)2  =  0.6;  (d)  (0\  =  0.3,  co2  =  0.7. 


to  one-third  of  that  of  the  original  structure.  However,  it  should 
be  noted  that  although  the  results  of  topology  optimization  pro¬ 
vide  good  guidance  for  better  material  distribution,  they  cannot  be 
regarded  as  fully  optimized  structures. 

5.  Discussions 

There  are  too  many  types  of  PEMFC  stack  products  to  list.  It  is 
impossible  to  completely  consider  all  the  variable  features  at  once 
owing  to  their  numerous  possible  combinations.  In  this  paper,  sev¬ 
eral  common  cases  are  considered  to  illustrate  the  influences  of 
design  parameters  on  the  optimization  results. 

5.1.  Effect  of  the  number  of  clamping  bolts 

Here,  two  additional  types  of  end  plate  structures  (Fig.  7),  con¬ 
taining  4  and  10  bolts  respectively,  are  studied.  Although  the 
optimized  structures  in  Fig.  7(a)  and  (b)  are  distinct,  there  are  still 
some  design  rules  available.  For  example,  to  keep  the  contact  pres¬ 
sure  as  uniform  as  possible,  the  end  plates  should  be  strengthened 
in  the  area  of  possible  large  deformation,  thus  increasing  the  local 
stiffness. 

5.2.  Effect  of  the  number  of  cells 

Taking  advantage  of  the  equivalent  stiffness  model  described 
in  the  literature  [12,13],  equivalent  stiffness  versus  the  number  of 
cells  can  be  rapidly  calculated.  According  to  Fig.  8(a)  and  (b),  the 
optimized  material  distributions  vary  with  the  number  of  cells.  This 
is  mainly  because  the  contact  response  changes  as  the  equivalent 
stiffness  of  the  stack  decreases  nonlinearly  with  an  increase  in  the 
number  of  cells  [13].  However,  the  optimized  material  distributions 
remain  similar. 

5.3.  Multi- objective  optimization  with  different  weighting 
coefficients 

The  two  weighting  coefficients  are  both  0.5  in  the  reference  opti¬ 
mization  model,  indicating  that  they  make  the  same  contribution  to 
the  final  objective.  In  a  practical  design,  the  two  coefficients  need 
to  be  adjusted  in  accordance  with  the  requirements  of  different 


problems.  For  example,  if  the  above  design  is  inclined  to  the  first 
single-objective  optimization,  the  first  corresponding  weighting 
coefficient  should  be  increased  appropriately.  Here,  topology  opti¬ 
mizations  for  another  four  groups  of  combinations  of  the  weighting 
coefficients  are  studied,  as  shown  in  Fig.  9.  Comparing  with  Fig.  5,  it 
is  seen  that  when  a  weighting  coefficient  is  increased,  the  obtained 
topology  of  the  structure  is  closer  to  its  corresponding  single¬ 
objective  optimization.  However,  it  cannot  be  simply  assumed  that 
a  similar  combination  of  the  coefficients  always  leads  to  similar 
optimized  topology  owing  to  the  nonlinearity  of  Eq.  (3).  Therefore, 
there  are  two  entirely  different  topologies  as  depicted  in  Fig.  9(b) 
and  (c).  In  most  circumstances,  several  trials  are  needed  for  engi¬ 
neers  to  obtain  the  best  combination  of  weighting  coefficients. 

6.  Conclusions 

In  this  paper,  a  multi-objective  topology  optimization  model 
of  the  end  plate  in  a  PEMFC  with  nonlinear  contact  boundary 
conditions  has  been  established  to  obtain  the  optimized  material 
distribution,  and  as  a  result,  maximize  the  stiffness  of  the  end  plate 
and  make  the  pressure  distribution  uniform  across  the  whole  cell. 
It  is  found  that  the  structure  designed  according  to  the  topology 
optimization  will  meet  both  the  requirements  of  low  weight  and 
manufacturability.  The  topology  optimization  is  the  basis  for  the 
shape  optimization  and  size  optimization  of  end  plates  and  could 
provide  good  guidance  for  designers  to  obtain  a  better  material 
distribution  at  the  early  design  stage. 
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